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1.  INTRODUCTION 


I 


The  computer  code  EMPFIT  was  written  to  give  a method  for  approximating  the  output 
waveform  of  an  electromagnetic  (EMP)  pulse  with  a relatively  simple  function  that  smooths  out 
small  numerical  variations,  but  is  easily  differentiated  and  Founer  transformed.  In  development, 
EMPFIT  was  coded  from  the  theoretical  work  of  Daniel  F.  Higgins.1  Also,  the  Fourier  transform 
code  used  in  EMPFIT  was  written  by  W.  Talmadge  Wyatt  of  the  Harry  Diamond  Laboratories 
(HDL).  As  can  be  seen  in  appendix  A,  EMPFIT  produces  excellent  results,  but  some  care  must 
be  exercised  when  choosing  particular  input  parameters  needed  to  run  the  code.  This  aspect  is 
dealt  with  fully  in  section  4. 

A quite  useful  code,  EMPFIT  requires  only  a minimal  number  of  input  data  points  to  find  a 
smooth  fit  to  a waveform  over  several  decades  in  time.  Normally,  digitization  of  a waveform 
requires  from  60  to  500  digitized  points  or  more,  depending  on  the  wanted  degree  of  accuracy 
and  refinement,  but  EMPFIT  allows  one  to  describe  the  waveform  with  as  little  as  15  to  30 
points.  This  aspect  makes  EMPFIT  most  attractive  where  the  user  desires  a smooth  curve 
through  data  points  and  a '‘nice”  Fourier  transform  (that  is,  low  noise,  smoothness).  Thus,  this 
advantage  makes  EMPFIT  a very  useful  tool  in  analyzing  and  describing  EMP  waveforms,  as 
well  as  other  data  that  require  smooth  curves  and  nice  Fourier  transforms. 

Regarding  the  operational  aspects  of  EMPFIT.  the  code  was  written  for  a Control  Data 
Corporation  (CDC)  6600  computer  system  at  the  Mobility  Equipment  Research  and  Develop- 
ment Command  (MERADCOM),  Fort  Belvoir,  VA,  using  standard  FORTRAN  IV.  The  code  is 
run  by  using  the  SCOPE  3.4.3  control  language.  A Houston  Instruments  Complot  Plotter  along 
with  the  appropriate  software  is  used  in  EMPFIT  to  obtain  plots  at  the  HDL  Woodbridge 
Research  Facility  (WRF).  Listings  of  EMPFIT  and  the  control  cards  necessary  to  run  the  code 
are  given  in  later  sections. 

In  the  following  sections,  the  general  theory  behind  EMPFIT  is  summarized  along  with 
some  problems  and  errors  encountered  while  coding  EMPFIT.  A detailed  discussion  of  the 
results  is  included,  as  well  as  numerous  examples  of  common  problems  and  solutions  to  the  most 
recurrent  problems.  The  output  options  of  EMPFIT  are  numerated,  a brief  description  of  the 
subroutines  of  the  code  is  given,  and  a detailed  section  on  preparing  data  cards  to  run  EMPFIT 
is  presented.  In  the  appendices,  a sample  run  and  examples  of  input  data  are  given. 

2.  CURVE  FITTING  THEORY 

A general,  short  description  of  the  theory1  for  EMPFIT  is  now  given.  Let  the  input  data 
points  (6,/)  for  1 s is  N be  given  that  describe  the  general  shape  of  any  waveform  /(r)  that  we 
wish  to  fit.  We  can  approximate  /(f)  for  / s f,  by  an  exponential  that  varies  as  e*.  and  for  t ^ ts 
we  can  approximate /(f)  by  an  exponential  that  varies  as  e~".  For  any  time  interval  t(sisi,  + | 
where  1 s i s N - 1 , we  fit  the  data  points  with  the  function 

~ 4)  + I ~ 0 
4+t  ~ 4 

■F  fii+i)(f  — 4)ft  — 4+i)  "F  Ci(f  — 4)(4+i  ~ /)*  + D,+I(fj+,  - t)(t  — l<)*. 


' Daniel  F.  Higgins.  A Method  for  Fitting  EMP  Waveforms  that  Facilitates  Calculation  of  the  Time  Derivative  and 
Founer  Transform,  Defense  Nuclear  Agency  Repon  DNA  .12} IT  (November  197}). 
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where 


//+,  ~f  f-fi-Ai 1_\ 

'*1+1  ~ l(  tf  ~ f|_i'  '/|+|  — tj 

-A/I"  + A/|*'  {^zl) 

(*<+i  — — *(-i) 


for  2 < / ^ )V  - l , 


for  2 s / s N - l , 


(2) 

(3) 


-A/I"  - A/|*'  P-^1) 

(*i  — ~ ti-i) 


for  2 s i rs  N - l . 


(4) 


A brief,  general  description  concerning  the  fit  of  equation  (1)  is  as  follows:  the  first  term  in 
equation  (I)  is  just  a linear  fit  between  the  rth  and  (i  + l)st  data  points;  the  second  term  is  the 
quadratic  correction  based  on  an  average  curvature  (/?,);  the  third  and  fourth  terms  of  equation 
(I)  are  used  to  insure  that  the  first  and  second  derivatives  are  continuous  at  the  data  points.  The 
values  A/|"  and  A/|*’  are  given  by 


A/|" 


_ ft+t  ft 

*.+  ■ - It 


_ ft  f- 1 

• I - 1,-x 


- 1,) 


(5) 


|(2?(  + B,. , )l,  /(_  i), 

A/|*'  = Bhi  - B 

The  exponential  functions  fitted  to  the  front  and  rear  of  the  data  points  (l,,f)  are 

f(t)  = A,e-  + Ase*“',  for  isih 
/(/)  = A,e~8'  + A4e~*fl',  for  t a ts, 

where 

A,  * (/,  - A,e*-')e-"'', 

A,  = 

[«*  (-  j ""*)  + 

*.  m~~~  + K*.  + *,)(*,  ~ t ,)  - «/„ 

h «i 

a*/,. 

A,  - [/  - A4e-***>*\ 

A4  , 

[fi  + P (-■  Y '')]<,~W“ 

~ j j + - <»-i)  — P/ni 

^4  “ 0%  ~ 4-  5j.|). 


(6) 

(7) 

(8) 

(9) 

(10) 

(ID 

(12) 

(13) 

(14) 

(15) 

(16) 
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The  special  parametric  values  for  a and  f3  are  discussed  in  section  4.  The  following  special  values 
for  fl, , BSs  C,,  and  DN  can  now  be  given  as 


*•  - *>■  [*■  * *•  (^A-*  (H-1) + *)] 


D"  6 (f*  - (20) 

Using  the  above  coefficients,  we  can  evaluate  /(/)  at  any  time  t.  The  function  given  in  equation 
(I)  is  continuous,  passes  through  the  data  points,  and  has  continuous  first  and  second 
derivatives. 

3.  PROBLEMS  AND  ERRORS  ENCOUNTERED  IN  CODING  EMPFIT 

In  the  process  of  coding  EMPFIT  from  the  theoretical  work,  several  typographical  errors 
were  found  in  Higgins'  report,1  which  are  enumerated  below. 

Higgins'  equation  (4)  should  read 


instead  of 


Higgins'  equation  (20)  should  read 


instead  of 


A,  « (/*  ~ A«e-“'’)e» 


A,  = (fs  - A4e-”'")e*\ 


Higgins’  equation  (2 1 a)  should  read 


instead  of 


^ = «[/(r)  + A,**-] 


dj/Xt)  m At)  + A,e~ 
dt  a a 


1 Daniel  F.  Higgins.  A Method  for  Fitting  EMP  Waveforms  that  FaciliUtes  Calculation  of  the  Time  Derivative  and 
Fourier  Transform.  Defense  Nuclear  Agency  Report  DNA  323IT  (November  1973). 
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Higgins’  equation  (21c)  should  read 


dm 

dt 


= -£[/</)  + iV-*®] 


instead  of 


d/\l)  = JM  _ A4e~tm 
dt 


B P ' 

Higgins'  equation  (34)  should  read 

A/.1*1  = -(ISC’.  + - O - (6C._,  + 18D„)(/„  - /„_,) 


instead  of 


bfn13'  = -(I8C.  + 6 £>„+,)('„+,  - tj  ~ (6 C.-,  + I8D,)(/,  - /„+1). 


Also,  the  following  corrections  were  supplied  by  Mission  Research  Corp. 
Higgins'  equation  (2)  should  read 


^ = / ■/»  _ fn  y»-l\ / I \ 

^ii+i  ~ fa-i'  'f«+i  ~ tm-v 


instead  of 

^ _ /■/»♦!  ~ ■/»  y»  ~ i \ 

f»+l  — 4i  f«  **  Ai-|'  'fa+1  ~ 

Higgins'  equation  (7)  should  read 


A yni  = 


fa+i  - /„  - f,-,  + ~ + Bn-,)(tK  - *„_,) 


instead  of 


4/i"  + ,lA  + - '■>  - M.  - 

Higgins'  equation  (12)  should  read 


C, 


**  - 3a*  [*,  + *,  {—■)/ («*  (V^)  + “)] 


6(/,  - /,)* 


instead  of 


C, 


*■ + 3a*  [*• + **  (^V1)/1  (“*  (S^) + a)] 


6 (/,  - /,)' 
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i 

i 


dn  = 


Higgins'  equation  ( 16)  should  read 

6UjV  — /jv— 1 )* 


instead  of 


Pv  = 


6(/iv  - /*_,)* 


There  was  also  encountered  a considerable  amount  of  difficulty  in  employing  the  Fourier 
transform  as  calculated  by  Higgins.1  This  trouble,  which  was  mentioned  by  Higgins,1  occurred 
when  computing  some  complex  exponential  terms  in  the  Fourier  transform.  As  it  turns  out,  there 
was  encountered  some  high-order  cancellation,  which  involved  the  exponentials  e1"'"  written  as 
(cos  wtm  + i sin  ojIJ  in  Higgins'  equation  (36).  This  round-off  problem  occurs  since  the  CDC  6600 
series  computer  has  only  14-digit  accuracy  in  single  precision,  and  information  in  the  sixth  and 
higher-order  terms  is  lost  when  the  cosine  and  sine  are  evaluated,  and  it  is  just  these  terms  that 
are  required  to  find  the  Fourier  transform.  Even  when  carrying  28-digit  accuracy  in  double- 
precision calculations,  the  high-order  cancellation  is  still  evident.  Although  Higgins'  dealt  with 
this  problem  his  report  did  not  provide  adequate  information  to  solve  this  difficulty.  Instead  of 
generating  the  necessary  coding  to  handle  the  problem,  the  employment  of  an  existing  Fourier 
transform  routine  was  decided  upon.  As  can  be  seen  by  the  various  plots  in  appendix  A,  this 
Fourier  transform  code  gives  very  nice  results. 


4.  DISCUSSION  OF  RESULTS 

Several  troublesome  difficulties  encountered  when  running  EMPFIT  must  be  overcome  to 
run  the  code  effectively.  First,  care  must  be  taken  when  choosing  the  data  points  to  describe  the 
waveform  of  interest.  It  has  been  noticed  that  the  spacing  of  the  points  is  somewhat  arbitrary, 
with  the  guideline  that  the  density  of  the  data  points  should  be  greatest  where  the  function  being 
fitted  varies  most  rapidly.  This  result  can  be  seen  in  figures  1 and  2.  Notice  how  figure  2 shows 
data  points  being  taken  where  the  slope  varies  the  greatest.  Figure  I is  an  example  of  poorly 
chosen  data  pairs,  in  that  not  enough  points  were  taken. 

Another  important  aspect  to  note  is  that  the  values  of  a and  p in  equations  (7)  and  (8)  are 
very  critical  in  the  goodness  of  fit  to  the  data  points  at  the  front  and  rear  of  the  waveform.  For 
waveforms  that  start  at  about  I shake  (1  shake  = 110'*  s)  and  end  at  about  1000  shakes,  some 
good  values  for  a and  p are 

« = 1.210*, 

P = 5.0104. 

These  values  are  variable  and  change  according  to  the  particular  data  points  of  interest,  but  these 
figures  have  been  found  to  be  the  best  in  fitting  the  exponential  functions  (7)  and  (8)  to  the  ends 
of  the  waveform.  It  is  important  that  the  last  few  data  points  used  to  describe  the  end  of  a 


' Daniel  F.  Higgins,  A Method  for  Fitting  EMP  Waveforms  that  Facilitates  Calculation  of  the  Time  Derivative  and 
Fourier  Transform.  Defense  Nuclear  Agency  Report  DNA  .X23 IT  (November  1973). 
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Figure  I . Curve  that  does  not  follow  data  points. 


waveform  should  be  decreasing  toward  zero  so  that  the  exponential  function  (8)  fits  a curve  to 
the  last  data  point  that  asymptotically  approaches  zero.  If  a value  for  a is  desired  other  than  the 
one  recommended  above,  then  a good  guideline  to  follow  is  to  choose  a such  that 


This  has  been  found  generally  to  give  reasonable  values  of  a and  facilitate  a good  fit  to  the  front 
of  the  waveform. 


A third  troublesome  point  that  happens  occasionally  is  the  fitting  of  the  peak  amplitude 
value.  It  sometimes  occurs  that  the  peak  amplitude  data  point  is  overshot  by  the  curve  being 
• fitted  to  the  data  pairs.  Then  it  has  been  found  that  the  peak  amplitude  point  is  matched  only  if 

the  data  points  are  chosen  very  judiciously.  A rule  of  thumb  to  alleviate  this  problem  is  to 
choose  the  closest  two  points  on  both  sides  of  the  peak  value  to  have  corresponding  time-change 
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Figure  2.  Solution  to  problem  where  curve  does  not  follow  data  points. 
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separations  and  approximately  the  same  amplitude  value.  This  correspondence  matches  the 
average-curvature  term  B,  on  the  opposite  sides  of  the  peak  amplitude  value,  except  that  the 
slopes  are  the  same  but  have  opposite  signs.  Hence,  the  curve  going  through  these  points  goes 
through  the  peak  amplitude  value,  and  the  slopes  on  either  side  of  the  peak  value  are  the  same, 
but  are  opposite  in  sign. 

One  final  remark  about  fitting  various  types  of  waveforms  concerns  fitting  a very  steeply 
rising  waveform.  A useful  procedure  is  to  choose  as  the  first  data  point  the  peak  amplitude 
value.  Then  with  a suitable  choice  of  a,  the  fitting  is  done  appropriately,  as  can  be  seen  in  figure 
9 (p.  15). 

These  four  problems  and  other  frequent  problems  are  summarized  in  table  I.  This  table 
outlines  specific  problems  and  gives  corresponding  solutions.  It  also  refers  specifically  to  figures 
I to  16,  which  show  the  problems  and  solutions.  Examples  of  correct  plots  from  EMPFIT  can  be 
seen  in  figures  A-l  and  A-2. 


II 


Figure  3.  Curve  that  does  not  follow  data  points. 


TABLE  I. 

Common  Problems  and  Solutions  When  Running  EMPFIT 

Problem  Solution  Figure 

Curve  does  not  follow  general  outline  Add  more  data  points  to  describe  trace  1,  2,  3,  4,  5, 
of  data  points  more  fully  6,  7,  8 

Waveform  is  steeply  rising  Pick  peak  amplitude  value  as  first  data  9 

point  and  then  choose  a accordingly 

End  of  waveform  is  not  smooth  and  Decrease  value  of  p one  order  of  mag-  10,  11,  12 
approaches  zero  with  too  great  a nitude  and  add  more  points  to  de- 
slope scribe  trace  more  ftilly 

Front  of  waveform  differs  markedly  in  Decrease  value  of  a one  order  of  mag-  13,  14 
sign  and  form  from  rest  of  waveform  nitude 
End  of  waveform  does  not  approach  Increase  maximum  time  to  be  plotted  IS,  16 
zero  (TMAX) 

Peak  amplitude  value  is  overshot  Pick  closest  points  on  either  side  of  — 

peak  amplitude  value  to  have  equal 
time-change  steps  and  approxi- 
mately same  amplitude  values 


i i-otxi  TuauMA-r 
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ELECTRIC  FIELD  STRENGTH- VERTICAL 


Figure  12.  Solution  to  problem  where  curve  does  not  approach  zero  nicely. 


ELECTRIC  FIELD  STRENGTH-VERT 1CRL 


MAGNETIC  induction-azimuthal 


Figure  16.  Solution  to  problem  where  end  of  waveform  does  not  approach  zero. 


5.  OUTPUT  OPTIONS 

The  following  options  are  available  when  running  EMPFIT.  The  specific  places  where  the 
card  entries  are  made  are  documented  in  section  7. 


a. 

b 


c. 

d 


Multiple  runs 

Plot  titles  and  axes  information  (supplied  in  program) 

(1)  T vs  Er — time  versus  electric  field  strength — radial 

(2)  T vs  Ev — time  versus  electric  field  strength — vertical 

(3)  T vs  B« — time  versus  magnetic  induction — azimuthal 

(4)  T vs  JR — time  versus  current  density — radial 

(5)  T vs  Jv — time  versus  current  density — vertical 

(6)  T vs  a — time  versus  air  conductivity 

Number  of  points  calculated  between  input  data  points 
Maximum  time  to  be  calculated 
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c. 

f. 


Fourier  transform — yes  or  no 

(1)  Starting  frequency 

(2)  Maximum  frequency  to  he  calculated 
Plots — yes  or  no 

(1)  Linear 

Jr°g  !°g  Available  on  curve  fit  plots  only;  Fourier  transform  plots  are 

(3)  Semilog  ' 

Abscissa  (log)  , dlways  log  log 
Ordinate  (linear)  J 

(4)  Minimum  value  of  ordinate  to  he  plotted  (does  not  plot  data  points  below  this 
value) — useful  when  plotting  linear  or  semilog  and  range  of  values  is  very  large 

(5)  Curve  fit  plots — three  plots  may  be  obtained: 

(a)  Data  points  (plotted  with  X's)  and  curve  through  data  points 

(b)  Curve  only 

(c)  Both  (a)  and  (b) 

(6)  Plot  and  axes  labels — may  be  entered  by  user  if  information  supplied  by  program  is 
not  sufficient 


6.  DESCRIPTION  OF  EMPFIT  AND  ITS  SUBROUTINES 

As  shown  in  appendix  A,  EMPFIT  consists  of  a main  program  and  16  subroutines.  A brief 
documentation  of  EMPFIT  and  its  subroutines  follows. 

EMPFIT  Main  program:  controls  the  running  of  the  program  and  determines  the  number 
of  runs  that  the  program  makes  by  reading  in  a run  parameter  card. 

JNPUTT  Subroutine:  reads  in  all  the  input  data,  except  plot  and  axes  labels,  if  desired  by 
the  user. 

CURFIT  Subroutine:  calculates  the  curve  fit  between  the  input  data  points. 

FORT  Subroutine:  calls  the  Fourier  transform  subroutine  FLINE,  determines  the 
delta  frequency,  and  separates  the  Fourier  transform  into  its  real  and  imaginary 
components. 

OTPUT  Subroutine:  performs  all  of  the  output  operations  and  also  calls  COLMNS; 

prints  out  the  input  data  points,  the  curve  fit  data  pairs  and,  if  wanted,  the 
Fourier  transform  data  points. 

COLMNS  Subroutine:  outputs  the  information  of  OTPUT  in  columns. 

PLOTT  Subroutine:  performs  all  of  the  plotting  calculations  and  calls  the  plotting 
package  ANAPAC  for  plotting  on  the  Houston  Instruments  Complot  Plotter  on 
the  Mohawk  Data  Systems  remote  batch  terminal. 

ANOTAT  Subroutine:  contains  all  the  plot  labels  and  titles  (sect.  5)  necessary  to  get  the 
curve  fit  and  Fourier  transform  plots. 

A1A3  Subroutine:  calculates  the  constants  A1  and  A3  of  equation  (7)  necessary  in 
fitting  the  exponential  to  the  front  of  the  input  data  points,  calculates  also  C,  of 
equation  ( 19). 

A2A4  Subroutine:  calculates  the  constants  A2  and  A4  of  equation  (8),  which  are  uted 
in  the  fitting  of  the  exponential  to  the  end  of  the  input  data  points;  calculates 
also  Ds  of  equation  (20). 
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FLINE 

BN 

CN 

DN 

DELFI 

DELF2 

ENTITL 


Subroutine:  calculates  the  Fourier  transform  of  the  curve  fit  data  points. 

Function  subroutine:  calculates  the  coefficients  Bt  of  equations  (2),  (17)  and 
(18). 

Function  subroutine:  calculates  the  coefficients  C,  of  equation  (3). 

Function  subroutine:  calculates  the  coefficients  D,  of  equation  (4). 

Function  subroutine:  calculates  equation  (5). 

Function  subroutine:  calculates  equation  (6). 

Subroutine:  reads  in  the  plot  label  and  axes  information  if  requested  by  the 
user. 


7.  DATA  INPUT  PREPARATION 

Input  data  cards  for  EMPFIT  are  prepared  in  the  following  manner.  Examples  of  input  card 
decks  appear  in  appendix  B. 

Column  Variable  Formal  Explanation 

Card  I : Multiple  run  card 

8-10  NRUN  13  Number  of  runs 


Card  2:  Plot  parameter  card 

10  IDENT  II 


20  . IFFT  II 

30  IPLOT  II 

40  ILNLOG  II 


Identifies  data  to  be  read  in 

IDENT  = I T vs  E„ 

= 2 T vs  Ev 
= 3 T vsBt 
= 4 T vs  J„ 

= 5 T vs  Jv 
= 6 T vs  <t 

Fourier  transform? 

IFFT  = 0 Yes 
= 1 No 

Plots? 

IPLOT  = 0 Yes 
= 1 No 

Plots  in  linear,  log-log,  or  semilog? 

ILNLOG  = 1 Linear 
= 2 Log-log 
= 3 Semilog 

Abscissa  (log) 

Ordinate  (linear) 

Note:  Fourier  transform  plots  are  always  log- 
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Column 

Variable 

Format 

Explanation 

41-50 

ORDM1N 

E10.3 

Minimum  value  of  ordinate  to  be  plotted  (all 
points  below  this  value  are  not  plotted  on  linear 
or  semilog  plots) 

60 

ICURV 

11 

Curve  fit  plots 

ICURV  = 0 One  plot  of  data  points  and  curve 
fit 

= 1 One  plot  of  curve  fit  only 

= 2 Two  plots — one  plot  of  data 
points  and  curve  fit,  one  plot  of 
curve  fit  only 

70 

Card  3:  Title  card 

IOT 

11 

Enter  own  plot  and  axes  labels? 

IOT  = 0 Use  plot  and  axes  labels  supplied  in 
program 

= 1 Enter  own  labels  on  cards  5,  6,  7 

1-80  TITLE 

Card  4:  Fitting  parameter  card 

8A10 

Title  or  subtitle 

1-7 

— 

— 

Blank 

8-10 

IPTS 

13 

Number  of  data  points  read  in 

11-17 

— 

— 

Blank 

18-20 

MPTS 

13 

Number  of  points  calculated  between  input  data 
points 

21-30 

TMAX 

E10.3 

Maximum  time  to  be  calculated  in  curve  fit 
calculations 

31-40 

ALPHA 

E10.3 

Used  to  fit  f[t)  = A,e“'  + A,e*"  to  front  of 
waveform  (good  starting  value:  a = 1.2E+8) 

41-50 

BETA 

E10.3 

Used  to  fit  J[t)  = Ate~m  + Ate  tB<  to  end  of 
waveform  (good  starting  value:  /3  = 5.0E+4) 

51-60 

OSTART 

E10.3 

Frequency  to  start  Fourier  transform  calcula- 
tions 

61-70  OMAX 

Card  5:  Abscissa  label  card 

E10.3 

Maximum  frequency  to  be  calculated 

Note:  Cards  5,  6,  and  7 are  used  only  if  IOT  = 1 
on  card  2.  If  IOT  = 0,  skip  to  card  8. 

1-10 

XTITLE 

AI0 

X label:  start  in  column  1 

22 


Column 


Variable 


Formal 


Explanation 


Card  6:  Ordinate  label  card 

1-20  YTITLE  2AI0  Y label;  start  in  column  1;  on  output,  ordinate 

label  is  only  in  A 10,  A2  format  instead  of  2AI0. 

Card  7:  Plot  label  card 


1-40  ATITLE 

4A10 

Plot  label;  start  in  column  1 

Card  8:  Data  card 

1-10  T 

E10.3 

Time  value  of  first  point 

11-20  F 

EI0.3 

Amplitude  of  first  point 

Card  9:  Data  card 

1-10  T 

EI0.3 

Time  value  of  second  point 

11-20  F 

EI0.3 

Amplitude  of  second  point 

Note:  Card  8,  card  9,  . . . are  repeated  with 
respect  to  the  number  of  input  points  indicated 
on  card  4,  [PTS. 

Note:  Card  2 to  card  8,  card  9,  . . . are  repeated 
according  to  the  number  of  times  identified  on 
CARD  1,  NRUN. 


8.  CONTROL  CARDS  FOR  EMPFIT 

The  following  SCOPE  3.4.3  control  cards  are  necessary  to  run  EMPFIT  on  the  CDC  6600  at 
MERADCOM.  All  the  underlined  entries  in  the  list  are  variable  and  must  be  entered  for  each 
user. 

EM 

TASK  (TNEM PW TRTS)  [User's  name] 

ATTACH,  AGO,  BINEMPFIT,  ID  = EM7I602. 

ATTACH.  LIBA,  ANAPAC,  ID  = EM71605.  MR  = I. 

LIBRARY  (LIBA) 

MAP  (PART) 

AGO. 

7/8/9 

[Data] 

Q/6/7/8/9 


23 


APPENDIX  A 

SAMPLE  RUN  AND  LISTING  OF  EMPFIT 


This  appendix  shows  a sample  run  of  EMPFIT  and  lists  its  main  program  and  subroutines. 
Figures  A- 1 and  A-2  are  samples  of  correct  plots. 


MAGNETIC  I N0UCTIC'<- AZIMUTHAL 


TJME  I SEC) 

Figure  A-l.  Sample  plot  from  EMPFIT. 
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P tOC  A AM  A«0  BLOCK  ASSICNN'NTS. 


• LOCK 

AOONISS 

If NATH 

Pill 

OAT!  PN0CSS4 

VON 

IfVfl 

MAKONANI 

CONNENTS 

/A/ 

111 

17 

INPUT 

130 

42  24 

ICO 

12/1 A/ 74  PTN 

4.4 

420 

AAAI 

1 

OPT-1 

/N/ 

ASIA 

A2S2 

/€/ 

1PA10 

A 

/O/ 

I2A1A 

SA7S 

/C/ 

/0SI1 

1 

/H/ 

TOSH 

1 

INPUT! 

20S13 

ISO 

ICO 

22/14/74  PIN 

4.4 

420 

AAAI 

UPT-l 

/f/ 

TO*A) 

10 

CVJNMT 

20A  IS 

2S4 

ICO 

12/1 A/7 A PTN 

4.4 

420 

AAAI 

OPT-* 

/ P ✓ 

21127 

SA  70 

inoi 

2101  7 

S2 

ICU 

I2/IA/7A  PTN 

4.4 

420 

AAAI 

OPT-1 

OTPUT 

27071 

A S | 

LCD 

I2/IA/7A  PTN 

4. A 

420 

AAAI 

OPT  • I 

COINNS 

27S42 

A 0 1 

ICO 

12/14/ 74  PTN 

4.4 

420 

AAAI 

OPT-1 

PLOT! 

301A3 

312 

ICO 

12/14/74  PTN 

4.4 

420 

4441 

OPT-1 

ANOTAT 

30A7S 

ISA 

ICO 

12/14/74  PTN 

4.4 

420 

AAAI 

OPT-1 

PlINP 

so  AS  3 

221 

ICO 

12/14/74  PTN 

4.4 

420 

AAAI 

OPT-1 

A1A3 

3107A 

111 

ICO 

12/14/74  PTN 

4.4 

420 

AAAI 

OPT-1 

a;aa 

3I20S 

IIS 

ICO 

12/14/74  PTN 

4.4 

420 

AAAI 

OPT-1 

IN 

31322 

AS 

ICO 

12/14/74  PTN 

A.4 

420 

AAAI 

OPT-l 

CN 

S140S 

AA 

ICO 

12/14/74  PTN 

4.4 

420 

4441 

OPT-1 

ON 

31  ASS 

AO 

ICO 

12/14/74  PTN 

4.4 

420 

AAAI 

OPT-l 

iCIPl 

mu 

AA 

ICO 

12/14/74  PTN 

4.4 

420 

AAAI 

OPT-1 

omi 

1IA0I 

30 

ICO 

12/14/74  PTN 

A. A 

420 

AAAI 

OPT-1 

/ 1/ 

IIA3I 

SA 

(NIITl 

1IAA7 

AA 

ICO 

12/14/74  PTN 

4.4 

420 

AAAI 

OPT-l 

PLOT 

sins 

7A2 

VI -l  ION 

10/tl/TS  PTN 

4.1 

PS93 

AAAI 

OPI-I 

NNOSK 

S27SS 

IS 

VI  -1 IBA 

OS/2  7/79  PTN 

4.2 

74  949 

AAAI 

OPT-l 

ONQP 

S27SA 

It 

Ul-llBA 

0S/1T/75  PTN 

4.2 

T43AS 

AAAI 

OPT-1 

NANKIN 

S2777 

AS 

Vl-lIBi 

0S/2S/7S  PTN 

4.2 

74  3 AS 

AAAI 

OPT-l 

local  s 

S30AA 

124 

Ul-llBA 

07/01/79  PTN 

4.3 

P393 

AAAI 

OPI-I 

INAIIS 

SA010 

SIS 

UL-llBA 

0B/1B/7S  PIN 

4.3 

P3V3 

AAAI 

OPT-1 

NOAH 

SA12S 

111 

VI -1  ION 

09/1  B/n  PTN 

4.S 

P393 

4441 

OPT-1 

• INN  A I 

3 AOS  A 

AS 

VI  -l  ISA 

•2/09/74  PTN 

4.3 

PS9S 

4441 

OPT-l 

DRAW 

3S1AS 

ISO 

Vl-llBA 

01/31/74  PTN 

4.3 

P393 

4441 

OPT-l 

CONPAIS 

SSS1S 

A 

Ul-llBA 

II /OS/79  PTN 

4.S 

P IV  3 

AAAI 

OPT-1 

NOOSK 

SSS2I 

20 

VI -l  IBA 

09/2  7/79  PTN 

4.2 

TVS  AS 

AAAI 

OPT-l 

IMONN 

SSSA1 

its 

VI  k IB  4 

09/04/79  PIN 

4.3 

'393 

AAAI 

OPT-1 

IIM 

mu 

140 

VI -A  IBA 

02/41/74  PTN 

4.3 

P393 

•441 

OPT-l 

• NANA 

SBB9A 

nn 

VI  -1104 

02/09/74  PTN 

4.3 

P34S 

4441 

OPI-I 

NUMf  A 

STMS 

ttt 

VL-1144 

•t/91/74  PM 

4.9 

P199 

4441 

OPT-1 

SYMBOL 

402  SO 

SAI 

VI -l  IBA 

4.4 

420 

AAAI 

OPT-l 

/ STP.fNO/ 

A0S7I 

1 
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29 

/«. 10./ 

AMIS 

ItB 

APPENDIX  A 


LOAD  OftR  (RMIT  CTlfA  LOAOft  1.1-420  12/14/74  11.10.24. 


• LOCK 

aooarss 

if  46TM 

r tie 

OATC 

AROCSSA 

TIB  IfVfl  HARDHAAE 

CONNCOTS 

• M1KT* 

407>O 

• 

u -r  oataab 

04/22/74 

tQORASS 

1.  1-420 

RCL  INIT ILAl I2AT10N  BOUT IBf . 

COBIO- 

40740 

44 

ii -root  a an 

•4/22/74 

COORASS 

1.  1-420 

COBOON  CPOEO  I/O  BOUT  IBIS  AND  CONSTANTS. 

facoof  • 

41014 

121 

sl-r  oataab 

04/22/74 

COBRAS S 

1.  1-420 

rOAOAT If  0 HAITI  IBIO  (OBI. 

f ««*«• 

4 1 1 S 7 

41 

H-roa7taa 

04/22/74 

CONRASS 

1.  1-420 

INITIAL 12f  CONSTANTS. 

t l»9UT. 

41220 

111 

si-rotTaaa 

04/22/74 

CONRASS 

1.  1-420 

CONNON  riOATINC  Ol/TRuT  CODE 

TOASTS- 

41411 

402 

SL-f  at  TRAM 

04/22/74 

LOOT ASS 

J.  1-420 

FOAIAAN  UiJECT  LlAAAAT  UTILITIES. 

incon- 

42111 

274 

Si-roaTtAa 

04/22/74 

COORASS 

1.  1-420 

COBOON  INRUT  FOANATT 106  COOf 

I arc  • 

42411 

140 

SV.-#0«TAAA 

•4/22/74 

COORASS 

1.  1-420 

roaOATTCO  A«A0  FORTRAN  aicoao. 

■ aora- 

41011 

444 

si-roa  taa« 

04/22/74 

COORASS 

1.  1-420 

OUT  RUT  ROB  NAT  IBTlARBHtn. 

ou?c- 

41447 

172 

Si  -7  OR  TRA4 

04/22/74 

COORASS 

1.  1-420 

fOANATTlO  HAITI  1 Cl Tt  AN  AUUAD. 

soar 

41441 

41 

Si  -r  OR  TRAM 

04/22/74 

COORASS 

1.  1-420 

CONRUTI  TMf  SOUAAI  tOP»  Of  >.  • Al  L . 

S TS • 1 ST 

41724 

42 

SL-rOATRAA 

04/22/74 

COORASS 

).  1-420 

NATH  LIIAART  UM  IL  f AAliA  ■ ‘ S * G • RAOCESSO 

ATOM 

44004 

10 

K-roaraAfl 

04/22/74 

COORASS 

1.  1-420 

RIAL  TO  INl(6fR  f IRONfAT 1 AT  ION . 

1 TOT* 

44014 

7 

Si-roa  TRAN 

04/22/74 

COORASS 

1.  1-420 

A 1 Al  TO  Al Al  f IRONINT IAT ILN. 

T LTIA- 

44024 

144 

Si-roa  taao 

04/22/74 

COORASS 

1.  1-420 

CONNON  FLOAT  INC  iNRUl  CONVlATft. 

f »T»»  ■ 

44201 

142 

Si-rOA  TRAN 

04/22/74 

COORASS 

1.  1-420 

CRACK  ARLISI  INC  fOANAT  f( A * r 3 1 • /K A AK | A . 

» oauTi . 

44441 

14 

Si-r  OR  T AAA 

04/22/74 

CONRASS 

1.  1-420 

K L NISC  . UTILITIIS. 

6ETf IT. 

44471 

42 

st-r  or  Iran 

04/22/74 

COORASS 

1.  1-420 

LOCATE  AN  1 1 T 6ITIA  A Mil  NANI. 

/ 10. BUT./ 

44411 

227 

iorr- 

4 4042 

114 

SL-roaTAAM 

04/22/74 

COORASS 

1.  1-420 

BINAAT  AIAO  fOATAAN  AICOAO. 

KUMa* 

44174 

404 

Si-rOATBAM 

04/22/74 

COORASS 

1.  1-420 

RAOCfSS  rOANATTEO  fOATAAN  IN’uT. 

2U  »2- 

44004 

201 

Si  -roa  taar 

04/22/74 

CONRASS 

3.  1-420 

BINAIV  HAITI  10ATAAN  A f C UA  0 . 

OUTCUN. 

44207 

144 

SL-f ORTRAN 

04/22/74 

CONRASS 

1.  1-420 

CONNUN  OJTRuT  CUOI 

achind- 

44141 

17 

SL-fORTRAM 

04/22/74 

COORASS 

J.  1-420 

raSlTlOA  R ILL  AT  1 1 1 1 NN 1 Nl -Of  - 1 Nf CAN  AT  ION . 

ClOCK- 

44422 

11 

Si-roa  TRAN 

04/22/74 

COORASS 

1.  1-420 

ACCESS  STSTIN  CLOCKS  FUA  fOATAAN. 

60  TOC  A- 

44441 

14 

Si-fOATRAA 

04/22/  74 

COORASS 

1.  S-420 

COORUTfO  60  TO  EA«i  RACCfSSOt. 

4 LUC 

44447 

71 

Si -ROB  TAAO 

04/2  2/74 

COORASS 

1.  1-420 

COORUTf  CONNON  AND  NA’JtAl  LOG AK 1 THNS . ORT- 

f l» 

44442 

74 

Si-r  OAT  AAA 

04/22/74 

COORASS 

1.  1-420 

f iRONf NT i al  Function,  f ti  ROnf*  *.  3Rt.au 

SUNOS. 

44447 

44 

Sl-r  ORTRAN 

04/22/74 

COORASS 

1.  1-420 

T A 1 6&NPN (TRIG  SINf  UK  CJSIM  '■  f A.  CRT-All. 

S TSAI  0* 

44  744 

1 

SL-fOATAAN 

04/22/74 

COORASS 

1.  1-420 

LIMA  AETHtlN  ST  S" A ID  AND  INITIALISATION  COO 

/coo. an/ 

44744 

4 

c to. aa 

44  744 

40 

SL-STSIO 

04/22/74 

COORASS 

1.  1-420 

zaor  .ao/ 

4 7014 

10 

apvf  .aa 

47024 

44 

4L-S7SI0 

04/22/74 

CONRASS 

1.  1-420 

an  .aa 

47110 

211 

Si -STS  10 

04/22/74 

CONRASS 

1.  1-420 

/ jaas.aa/ 

47141 

11 

/acac  .an/ 

4 7144 

1 

/oars. co/ 

47147 

1 

zoara.ro/ 

47140 

7 

oar  a.aa 

47147 

214 

SL-STSIO 

04/22/74 

CONRASS 

1.  1-420 

/ TCR4.au/ 

4 7424 

1 

zauT.ro/ 

4 7424 

. 1 

aur.sa 

47444 

1142 

SL-STS  10 

04/22/74 

COORASS 

ft.  *-420 

haa.so 

41214 

240 

S4.-JTS  10 

04/22/74 

COORASS 

J.  1-420 

/asr.ro/ 

41474 

7 

cisr.aa 

41404 

24 

SL-STSIO 

04/2  2/74 

COORASS 

ft.  1-420 

/ 6f  T.tT/ 

41410 

4 

•TftT.lt 

41445 

Ilk 

jL-iuia 

04/22/74 

COlf  ft  15 

J«  S-420 

NfOl.lt 

41441 

142 

SL-STSIO 

04/22/74 

COBRAS  S 

ft.  1-420 

zsarL.ro/ 

42011 

7 

safL.ia 

92022 

47 

11-STilA 

•4/22/70 

CJBR4S5 

J.  1-420 

c aa.aa 

42c  7 1 

404 

SL-STSIO 

04/22/74 

COORASS 

1.  1-420 

c*ara.so 

42474 

7 

SL-STSIO 

•4/22/74 

COORASS 

1.  1-420 

jsua.aa 

42404 

44 

ii  -snip 

04/22224 

COBRftll 

A.  S-420 

oafo.so 

41471 

242 

SL-STSIO 

04/22/74 

COBRASS 

ft.  1-420 

LOAO  OAR  - 

roar  it 

CV9EA  LOAOtA  I 

BRCI.M 

51051 

14 

SL-SVSIO 

04/2  2/74 

COORASS 

1. 

1-420 

/RUT. AT/ 

5104? 

11 

ALfA.RN 

51100 

42 

SL-STSIO 

04/22/74 

COBRASS 

ft. 

1-420 

cisr.so 

11142 

152 

Si -STS 10 

04/2  2/74 

COORASS 

ft. 

1-420 

/CLSV.RO/ 

91274 

7 

C 1ST. SO 

51101 

121 

SL-SVSIO 

04/22/74 

COBRASS 

ft. 

1-420 

ZAIH.ro/ 

51N2A 

7 

BIN. SO 

51415 

SI 

SL-STS 10 

04/22/74 

COBRASS 

J. 

1-420 

zcfT.ro/ 

51444 

? 

Z6fT.iT/ 

51475 

II 

6FT.S0 

51504 

1015 

SL-SVSIO 

04/22/74 

COORASS 

ft. 

1-420 

I .so 

54541 

101 

U-SVSIO 

04/12/74 

CONRASS 

1. 

1-420 

H .SO 

54444 

50 

SL-SVSIO 

04/2  2/74 

CONRASS 

1. 

1-420 

r su.sa 

54714 

104 

SL-STS  lb 

04/22/74 

CONRASS 

S. 

1-420 

SRI  .to 

95022 

17 

U-OUCLfUS 

19/12/74 

COORASS 

1. 

1-420 
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APPENDIX  A 

• cast  an*  tip  mvitDNNiNT 

INPUT  DATA  • HAUHIIL  INDUC  I ION-42  INulMAA. 


pna*  1.2001*0*  im* 

N»M  Of  INPUT  DATA  POINTS* 
TIN!  ISICI  *-AI INN THAI 

* .000  f *0A  MAX IN UN  I INI  TO  01 

JA  NUMt*  Of  POINTS  CilCUtATfO 

TINA  < SCC 1 a-AIiaUTMAl 

CA1CUA.ATI0*  2.0001  -OA 

BITWIIN  INPUT  DATA  POINTS*  A 

7 INI  ISCCI  i-AI INUTNAl 

1 .71*1-07 

-*.a«ii-oa 

A.iaof  -07 

-O.A*  11-01 

A . 3001 -OA 

-A.  10 Of -01 

1 .AAAI -07 

-1  .*661  -06 

*.210f  -07 

-7.*A  71-01 

*.2*71-0* 

-3. *0*1  -01 

2.  077l>07 

-6. *6*1-03 

7 . 1A7(  -07 

-*•**31-01 

A. 3*11-06 

-2.3911-01 

2.27*1-0? 

-1.2611-42 

B.  1121-07 

■*•1311-01 

1.2111-0* 

-*.0311-02 

2.  AMI-07 

•2. *171-02 

O.OOAI  -07 

-*.l**f-01 

I .32*1-0* 

-A  .*0*1-02 

2.7*91-0? 

-9.29*1-42 

1 .00*1  -0* 

-*.2  711-01 

l.ASAf-O* 

-*.6*11-02 

2. tail -07 

-7.0171-02 

l.lOAf  -0* 

-*.10*1-01 

l.SOAl-O* 

-*.**21  -02 

3.0211-07 

-1.2701-01 

1 .7*11  -0* 

*.23*1-01 

2.7731-0* 

-6.1*71-02 

7.  3191-0? 

-7.7101-01 

2. *201 -OA 

-*  • A 1 11-01 

3. **61-0* 

-A. 7661-02 

9.0*11 -07 

-1  .AMI  *00 

l.OAlf-OA 

-4.2001-01 

* .21*1-0* 

-A  ,A07£ -02 

9. 0*71*07 

-1.0*71*00 

9.45*1  -OA 

-A.OAlf-O  1 

7.431! -0* 

-*.10*1-02 

A. 1001-07 

-1 .0001*00 

Auoaia  Of 

POINTS  • 1A 

PIT  (AKUiAT  IONS 

AIAJA  VAlUl 

Of  010INATI  TO  11 

PlOTTfD*  -1.0001*09 

T|N|  (SKI 

1-Af INuTHAi 

UNt  ISCCI 

• -A2INUTMAI 

TIN!  ISICI 

•-A1INUTHAI 

1 .37*1-07 

0. 

6.0771-07 

-1.0091*00 

6.1021 -0* 

-6.5791-01 

1 .6*31-07 

-3.9AM-0* 

6.1001  -07 

-1.0001*00 

6.2*01-0* 

-6.6621-01 

1 . *121 -07 

-1.1*31-0* 

6.1*61  -07 

-9.92*1-01 

6 . 1641 -06 

-6.3001 -01 

1 .*111-07 

-1.91*1 -0* 

6.2121  -07 

-9.9911-01 

6.5ITI-06 

-*. 12*1-01 

1.6*91-07 

-3.7911-06 

6 . 2**1  -07 

-9.101-01 

6.7961-0* 

-3.6671-01 

1.7111-07 

-9.6111-06 

6.3261  -07 

-9.52*1-01 

6.9)21 -06 

-3.7671  -01 

1 .7*21-07 

-*.01*1-06 

6.3101  -07 

-9.6931-01 

9.1091-0* 

-3.5*51-01 

1 . 7851-0? 

-3.6311 -0* 

6.  **71  -07 

-9.1)31-01 

5.7171  -06 

- 3.6091 -01 

I .*  191 -07 

-1.9*11-0* 

6. 7261  -07 

-9.7)91-01 

5.5011-06 

-9.1621-01 

1.1*21-07 

-S. 1791-05 

6**991 -07 

-*.1191-01 

5.7191 -06 

-2.6861  -01 

1.11*1-07 

-1. *661-06 

9.0*71-07 

-7.9261-01 

5.9)01-06 

-2.79  11-0  1 

1 .9231-07 

-*.7321-06 

*.2391  -07 

-7.99  71-01 

6.1661-06 

-2.9891-01 

1.9*01-07 

-1.9021-03 

9. *2  31-07 

-0.9121-01 

*.35*1-06 

-2.96*1-01 

I .9911 -07 

-2.9321-03 

*.0101  -07 

-0.61*1-01 

7.90*1-06 

-l .5951-01 

7.03*1-07 

-3. *711-03 

6.39*1  -07 

-*.03*1-01 

1.6591 -06 

-1 .0*61 -01 

2.0721-07 

-6.I6A&-QJ 

4. 7C1L.-02 

-1*J9  61-01 

9.1091-06 

-7. *7*1-02 

2.1131-07 

-*.1*11-03 

7.1671  -07 

-4.9691-01 

1.09*1-09 

-♦.2361-02 

2.1*61-07 

-7.90/1-0) 

7. SAW  -07 

-4.6701-01 

1.2111-05 

-5  .01*1  -02 

2.  19*1  -07 

-9.91*1-0) 

I.  *331-07 

-4.9191-01 

1.23*1-09 

-6.90)1-02 

2.23*1-07 

-1 .0*21  -02 

7.76*1-07 

-9.3771-01 

1.25*1 -05 

-6.7901 -02 

2.27*1-0? 

-1.2611-02 

7.9391  -07 

-9.9*21-01 

1.2AII-09 

-6.70A1-02 

2.3201-07 

-1 .6*61-02 

l.LMt  01 

-3.JJ1L-01 

1 .J05L-05 

-6. *671-02 

2.9*91-07 

-1. *911-02 

• . 3294  -07 

-4.9921-01 

1.32*1-05 

-6.6091-02 

2.6091-07 

-1.9*11-02 

1. Sill  -07 

-9.9931-01 

1.99*1-05 

-6.9*11-02 

2 •***!  -07 

-2.2*71-02 

1.  TIM -07 

-9.99  71-41 

1.3791-05 

-6.52H-02 

2. *911-07 

-2. *171  -02 

1.90  31  -07 

-9.9*91-01 

1 .6091 -09 

-6.6601-0  ; 

2. *6  71-07 

-2.9921  -02 

9.09*1  -07 

-*.96*1-01 

1.6901-09 

-6  .6661-02 

>.*971-07 

9. 21 M -02 

-5,3) U 'fli 

1.69*1 -05 

-6.65*1-02 

7. AAAI -07 

-3.9771-02 

*.*•21-07 

-9.32  91-01 

1.6*61-05 

-6.65)1-02 

7. *9*1-07 

-6.9191-02 

9.67*1-07 

-9.90  7I-OI 

1.9121-09 

-6.6591-02 

2.76*1-07 

-9.29*1-02 

9.1*71-07 

-9.2901-01 

1.5601-05 

-6.6601  -0  2 

2.7731-07 

-9.72*1-02 

1.00*1-0* 

-*.2731-01 

1.56*1-05 

-6.6701-02 

7. *001 -07 

-*.2191 -02 

1.02*1  -O* 

-9.2*  71-01 

1.59*1-05 

-6.6*21-02 

7.12*1-07 

-*. 711102 

1.0*91  -06 

-9.2*  11-01 

1.J1U-05 

-6.51*1-02 

2. ***1-07 

-7.3*11-02 

1.0*91-0* 

-9.2291-01 

2. 0*71-49 

-4.6791-02 

2 . *• 11 -07 

-7.9)71-02 

1.0*61-0* 

-9.1101-01 

2.9021-49 

-4.7*21-02 

2.9111-07 

-*.3*31-02 

1.10*1-0* 

-9.19*1-01 

2.53*1-09 

-4. *611-02 

2.93*1-07 

-*.  *601-02 

1.2391 -0* 

-9. 1601-01 

2.7711-05 

-6.8971-02 

7.96*1-07 

-9.9291 -02 

1 . 3*31  -0* 

-9.13)1-01 

2.950i  -05 

-6 .6A9f -07 

7.9931-07 

-1.0  731  -01 

1.6921  -0* 

-9.19*1-01 

1. 12*1  -05 

-6.1691-02 

3.0211-07 

-1.2791-01 

1.  *221-0* 

-9.1991-01 

3.3091-05 

-*.•)•!  -02 

3.0*11-07 

-2.0911-01 

1. 7911-0* 

-9.2991-01 

3.6791-49 

-*.*021-02 

3.1*01-07 

-3.2711-01 

1.90  71  -0* 

-*.2921-01 

3.69*1 -45 

-6.76*1-07 

1.7001-07 

-A.AAOf  -01 

2.0*21  -O* 

-9.92*1-01 

1.9921 -05 

-6.70*1 -0. 

3.7*91-07 

-*.19*1-01 

2.21*1  -0* 

-9.1*  71-01 

6. 1071  -0* 

-6 .6  181  -0. 

3.3191-07 

-7.71*1-01 

2.37)1  -0* 

-9.3991  01 

6 . * 1)1 -05 

-6.5661-02 

3.  !**f  -07 

-9.9061-01 

2. *29!  -0* 

-9.61  11-01 

6.65*1-05 

-6.5031-07 

3.6*31-07 

-1 .1791*00 

2.6311  -0* 

-9.60*1-01 

5 .2161 -05 

-6.6971-02 

1.9191-07 

-I .29)1 *00 

2.7961  0* 

-9.99  91-01 

9.7991-05 

-6.6061  -02 

3 . **Af -07 

-1  .62  31  *00 

7. *3*1  -06 

-9,96*1-01 

6.72*1 -US 

-6.7261 -02 

7.6*11-07 

-1 .6**1 *00 

2.9)91  -06 

-9.93*1-01 

6.6661  -05 

-6.7121  -02 

3.7201-07 

-1.6*61 *00 

) .0611  -0* 

-9.2991  01 

i.un-05 

6.5061 -02 

3.7*71-07 

-1.3771*00 

3. 1*61 -0* 

-9.26*1-01 

7.61*1-09 

-6.1091-02 

1.0*11-07 

-1.7*91 *00 

1.2*71-0* 

•*.1*61-01 

1. 0111 -06 

1 .629*  -02 

3.9201  -0  7 

-1.11*1*00 

».*10l  06 

-9.11  31-01 

1.25*1-06 

-9.07*1 -01 

1.9171-07 

-I  .0*71  *00 

1.  Mil  -06 

-9.0)61-01 

1 . S061  06 

•1.6061 -0. 

A. 0101-07 

-1.0601 *00 

1.6161 -06 

-4.96*1-01 

1 . 79  11-06 

-6. 6001  -06 

6. f 121  -Of 

-1.0211*00 

i.*0*t  06 

-4. *361-01 

2 . JOwt -06 

-1 .2*31*06 

6.0*91-07 

-1.01*1*00 

3.99  91-0* 

-6.7101-01 

2.26  71  -06 

-9.7211-09 
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APPENDIX  A 


»OjAI(R  IMISMI  CAUUlATKmS 

S’  UMtl  6RI  00l  NC  Y I Hit  1 1 1 ■ 1.0001*04  Dtllft  NlOUfXT*  1.0*71*00 

mtifjrn  mt^LlSd  TO  »r  (AllOlATfD*  ).000(*07 


• ruvimy  anhituoi  tftfttcacr  muruoi 

• 141  IMUNUT  REAL  IMAGINARY 


I .0001 *04 

i .4/3E-04 

1.0747-04 

). 7)71*0) 

1.3017-0  7 

1.  0107-07 

1 .04/1  *04 

-7.  7)71  -04 

1 .0707-04 

).)))f*0) 

1.4)77-07 

1.70)7-07 

; .004»*04 

-i.6777  -04 

1 . 006 ( - 04 

4.2741*03 

1.49)1-07 

1 . )04t -07 

1 . 1441 *04 

“7  • • 7)  E -04 

1.17)7-04 

4 .)  70  7 *0) 

1 .>>43-07 

1.4477-07 

J .7007*06 

-7  • )70 f -04 

1.1761-04 

4.6747*0) 

1.  >973-07 

1 .3013-07 

1 . 7441 *04 

-7.447E  -04 

1 .7497-04 

7.194f*0) 

1.61)1-07 

1.16)7-07 

I.IMC«|4 

-7.444E -04 

l .3137-04 

7 .>91  7 *0) 

1 .6)71-07 

1.0)07-07 

1 . 1 741 *04 

-7.47)1-04 

1.3647-04 

7.6677*0) 

J .6997-07 

9.1)17-06 

1 .4197  *04 

-7  .RRIE-04 

1 *4  )0  E -04 

6.7467*0) 

1.7101-07 

7.2447-06 

1 .4041*04 

-7.8  3)E  -04 

1 .)077-04 

9 .4)71*0) 

1 .6097-0? 

3.9207-06 

1 .4  741*04 

-7.7711-04 

l.))77-04 

9.0)47*0) 

1 .7773-07 

4. 6177-08 

1 .4441 *0* 

-7 . 70  I E -04 

l.)67t-04 

). 4)4  1*0) 

1.7711-07 

i .6467-06 

1 . 726f  *04 

-7.4317-04 

1 . ) 4 77 -04 

9.6947*0) 

1 .0777-07 

1 . 1097-00 

1 . • 041 *04 

-7.)4)E-0* 

l .))) 1 -04 

1 .03)7*04 

1.6497-07 

-4.7211-09 

1 . #4  IE  *04 

-7.)0)E-04 

l.)41|-04 

l .01)1*06 

l.»)i-07 

-2 . 164 1 - 08 

i .9797*04 

-7.44)f -04 

l.)007-04 

1.1)47*04 

1.3  701-0  7 

-3.7707-06 

7.0711*04 

-7.4)41-04 

l.)467-06 

1.1077*04 

1.6377-07 

-). 3337-04 

7.  1471*04 

-7.41)7-04 

1. >44  0-04 

1.2421*04 

1.3271-07 

-A.  0451-06 

7.7467 *04 

- 7 . 4041  -04 

1.4777-04 

l .300  7*04 

1.2011-07 

-6.76)7 -Of 

7.47)1 *04 

-7.3)41-04 

1.4477-04 

1 .3407*04 

1 .0323-07 

-9.3767-01 

>.4*4»  *04 

-7  . )•  73  -04 

1 .4)43-04 

1 .4747*04 

• .6407-0* 

-I  .06 9»-07 

7.3997  *0* 

-7.34)7  -04 

1 .4447-04 

1 .4901*04 

6.94)7-08 

-1.1)97-07 

7. 7707  *0* 

-7.3771  -04 

1.7)67-04 

1 .3)97*04 

4.6767-06 

-1.2147-07 

7.1471*04 

-7.7A7E  -04 

1 .0717-04 

1 .4)77*06 

2.7)17-04 

-1.2447-07 

7. 4*01 *04 

-7.147E-04 

1 .0747-04 

1 .706  7*06 

3.70)7-09 

-1.2207-07 

).i  m *04 

-7.1071  -04 

1.41)7-04 

l .70  77*06 

-1  .4)37-00 

-1 .1 767 -07 

). 74*1*04 

-7.07)1 -OA 

1.4)07-04 

1 Ulilt  04 

-3.6007-00 

-1. 0137.-07 

3.41)7 *04 

-1.4)11-04 

1.4047-04 

1 .9)07*06 

-).6413-00 

-6.4961-06 

3 • ) 747 ♦ 04 

-1  .R44E-04 

7.07)7-04 

7.0697*06 

-6.9303-06 

-7.3697-08 

3.7411 *04 

-1.747E  -04 

7.0)47-04 

7.1467*06 

-0.0427-06 

i.  9727  -06 

3.91)7 *04 

-1 .4)4#  -04 

7. 044 f- 04 

7.744  7*04 

-f. 7227 -Of 

-3.9JD7 -C* 

4.0)71*04 

-1  .)  77|  -04 

7. 046 7-04 

7 .3467*06 

-•.9)67-01 

-1.7697-08 

4.7677*04 

-1.3307 -04 

7. 0411-04 

2.6)07*06 

-■  .6  7i£  -0  6 

3.0647  -09 

4.4R7»*04 

-1.3)41-04 

7.1747-04 

2. >777*06 

-7.9673-00 

7.2277-06 

4.4041*04 

-1.7771-04 

7.13)7-04 

7.697  7*06 

-6.6093-00 

3.9)27-06 

4.9147 *04 

- 1 . 1)41  -04 

7.1747-04 

7.0177*06 

->.3073-00 

3.3147 -06 

).  |4)f  *04 

-1  .0441  -04 

7.1UE-C4 

2 .946  7*06 

-3.3767-06 

0.2661 -06 

4.  )R>f  *04 

-4.374E  -07 

7.040 f-04 

3.06)7*06 

-1  .3)67-06 

6.71)7-06 

4.4IIE  *04 

-1.7171-07 

7.07)1-04 

3.7797*06 

4.7817-09 

6.3927-08 

4.0)1  *04 

-7.0401-07 

7 .0  34 1 -04 

3.3797*06 

2.3977-06 

3.66)7-01 

4.1 70# *04 

-4 .004# -0  7 

1 .4947-04 

).5)4r*06 

4.0047-06 

4.3817 -06 

4.4)71*04 

-4. 447#  -07 

1 .9497-04 

3.7013*06 

>.1137-06 

2. 7967-06 

6.7)77*04 

-3.4477-07 

I .69/7-04 

).67)r*06 

>.31)7-06 

7.4037-09 

/.r  7«i  *04 

-7.47)1-07 

1 . 6 37  £ -04 

4.0)41*06 

3.1407-08 

-1 . »09f-08 

7 . 40  It *04 

-7. 0431  -C  7 

1 .7417-04 

4.7477*06 

4.0247-08 

-).042t-0S 

7.  74)f  *04 

-1.1)3#  -07 

1 .4697-04 

4.4407*06 

7.3791-08 

-4.2047-08 

R. 1041 *04 

-3.64)7-08 

1 »))6  7 -04 

4 .646  3*04 

1.3337-09 

-4.3697  -08 

• .44)7  *04 

7.)0RE -04 

l .)  1 ■ 7 - 04 

4.6637*06 

-1  .6267-08 

-4.06)7-0# 

• .6787  *04 

• .44)1 -0» 

1 .4747-04 

3 .0691*04 

-3.1 777-04 

-7 . 74 -f-u* 

4.7)11*04 

1 .3477  -07 

1.3347-04 

3. 3701*06 

-3.9767-08 

-4.291* - J9 

).  77JE  *04 

1.  7)67-07 

1.7477-04 

). >741*06 

-3.0407 -08 

1.219* -07 

1 .0167*0) 

7 .0))# -07 

1.1)47-04 

3 .6)31*06 

-2.1)07-08 

2 .8267 -08 

1 .04M*C) 

7.7/77-07 

1.0777-04 

6.1047*06 

-3.4)47-09 

1.4947 -09 

1.11)1*0) 

7 .4331  -0  7 

). 64)7-07 

6.3093*06 

I .6113-06 

2.9497 -Of 

1.1441 *0) 

7.)0)#-07 

4.1007-07 

6.6661*06 

7.91)7-08 

1.  *2)7-06 

1 *771#  *0) 

7 . )1  77-07 

6.3)47-07 

6.9971*00 

7 .9337  -07 

- 7 .46)7 -3? 

1 .7  771  *0) 

7.4)47-01 

7.47)1-07 

7.37  3 7 *06 

1.6697-06 

-2.1 707 -Of 

1 * ) ) 7f  *0) 

7.347E-07 

6.9767-07 

7. 664 r *06 

-3.6807-09 

-2.7127-06 

1 . )))#  *0) 

7.7441  -07 

4.3717-07 

6.070  7*04 

-2.0317-0# 

-1.3  7)7-06 

1 .4441*0) 

7.0441-07 

). 64)7-07 

4 .39)7*06 

-2.3847  -06 

1.6017-09 

1 .))7»  *0) 

1 .4)71-07 

). 3431-07 

3 .7441  -04 

-1.1*07-08 

».92)e-)6 

1 .404*  *0) 

1 .4477  -07 

4.6477-07 

9.19)7*06 

7./ 167 -09 

1 .9)*t -0* 

1.4/91*0) 

1.4401-07 

4.4747-07 

9.6711*06 

1 .6467-06 

3 . 06  It -09 

1 . 7)  7f  *0) 

1 .7417-07 

4.3)77-07 

1 .0073*07 

1 .7967-06 

-1  .1)67 -06 

1 .«»•» *0) 

1 .04)7-07 

4.13)7-07 

1 .0)47*07 

-3.1017-09 

-1 .3267 -06 

1.4741*0) 

4.77)1-04 

3.9777-07 

1 .10)7*07 

-1.3)17-06 

-i.7171 -09 

7.01  IE  *0) 

7.4707-04 

3.6447-07 

1 .1)47*07 

-7.6697 -09 

9.  * 1 ) 7 • )• 

7.I07E  *0) 

4. )4  7| -06 

3.73)7-07 

1 .706  7*07 

3.1677-09 

9.0077  -,)V 

7 . 70)t  « 0) 

).  44)1-00 

3.4147-07 

1.2061*0  7 

8.3497-09 

-1.6691-39 

7.30RE *0) 

). 4)47-04 

>•>1)7-07 

1 .37)7*07 

•2.3793-1 1 

-7.2707-09 

7.41)1*0) 

4.4007-04 

1.4747-07 

1 .3047*07 

-3.6193-09 

-6.791! -10 

7 . )7»t *0) 

4.44)1-04 

3.3477-07 

1 .4697*07 

-9.7)03-10 

4.3677  -09 

7 . 44)E  *0) 

4.7477-00 

3.7  797-07 

1. >167*07 

3 .3693  -09 

) .99  »<  -1C 

7.  74iE  *0) 

4.7741  -04 

J.7  347-07 

1 . >077*07 

-6.4)07-1) 

-2.7)11-09 

7 . •) 71 *0) 

4.4711-04 

3.10)7-07 

1 .6011*07 

-2.0207-09 

7.794i-ll 

3.0)71*0) 

). 0477-06 

3.1777-07 

1 .7)07*07 

1.42)7-09 

1.2)37  -09 

).l 7)1*0) 

). 4447-06 

3.0)07-07 

1 .019  7*0  7 

3.1077-10 

-1 .6)07 -09 

1.3711*0) 

4.1177-00 

7.9477-07 

1 .9047*07 

-1  .0713-09 

4.7)27-10 

3.4  7)7  *0) 

4 . )4  7#  -Of 

7.6777-07 

1 .9971*0? 

1 .4663-09 

1 ..947 -09 

3.4371  *0) 

7.044*  -06 

7.6177-07 

7.00)7*0  7 

3.3497-10 

-1.40)7 -0« 

3 • *041 *0) 

7.7711-00 

2.7)91-07 

2.1071*0? 

-1.6)17-09 

7.) 771-10 

). 9647*0) 

6.4477-06 

7.4493-07 

7.70)7*07 

1.40)7-09 

7. 4)ff -|0 

4.14)1*0) 

4.0607  -«• 

7. >>43-0? 

7.3907*07 

-6.7793-10 

-1.3661 -04 

4.»*>(*0) 

4.40)1-01 

7.4771-0? 

7 .)0l  3*0  7 

-3.7097-10 

1.9041-09 

4 . ) 44  E *0) 

1.0447  -0  7 

3 .3447-07 

7 .4177*07 

6.6117-10 

- 3 . 1 )) f -10 

4. 77RE  *0) 

1.1)17-07 

7. 7477  -07 

7.7397*07 

-6.7627-10 

-3.4647-11 

4.4T)lt*0) 

1 .'1M-0T 

M»FW 

?.# 4tl  + r 

I.IIOf-lO 

*.4?W-!1 

4.7341*0) 

1 .7071-07 

7 .0)03-07 

3.0003*07 

1.2003-10 

->.3607-11 

APPENDIX  B 
SAMPLE  INPUT  DATA 


This  appendix  shows  examples  of  input  card  decks  of  EMPFIT. 


1 

2 5-1.0007*05  2 

(NVlft3M*!4T  III 

49  4 1.0001-02  1.2001*09  5. 000140?  l.000t*04  1.0001*07 

I .26)7-07-6.1617*00 
1 .4 147-07-1. 724f*01 
1 .S44f-07-4. #157*01 

1 .7721-07-1  .44*00? 

I .976I-07-*. 7741*02 
2#l  •1007-7. #7500? 

2.4657-07-1.6267 *01 

2 . 550007-2.  5547*01 
2 .4507 “0  7-  5. 9571*01 
1.0197-07-2.440*04 
1 .1171-07-1.1147*04 
4.0#07-C7-?. 7427*01 
4.4917-07  7.0507*0? 

5.4797-07  1.2417*01 
4.4417-07  4.6927*01 
7 .1077-07  4 . #097*01 
I .0647-0?  4.70#i*0l 
4 .4747-07  4.7477*01 
9.7407-07  4.5247*01 

1.1417- 06  1.4627*01 
1 .5647-04  1 .4457*01 
1.7517-06  1.7417*0? 

1 .4907-06-4.9297*02 
2.0C07-06-2  . 2707*01 

2 .5291-04-4.0207*01 
2.4C07-C6 -6. 7167*0: 

2 .4007-06-7.44  77  *01 
4 .5047-06-1  . 0227*04 
5.7957-06-1 .14?{*04 
1 .1057-05-1.1207*04 
1 .2117-05-1  .IM!*04 
1 .1247-05-1  .1247*04 
1 .4567-05-1.1417*04 

1 .7517-05-1  . 1 747  *04 

2 .5107-05-1  .2977*04 
1.1147-C5-I .1727*04 
5 .7447 -05  -1  .4777*04 

9. 1417- 05-1.  5017*04 
1 .0007-04-1  .4997*04 

1 .100I-C4-1  .4707*04 
2.2007-C4-l.a9'#7*04 
2.5507-04-1  .0451*04 

2 .4501-04-1 .0007*04 
4.200f-C4-5.?00»  *01 
5.0007-04-4.4007*01 
6.5007-04-1. 4007*01 
9 .0007 -06- 2.0007*01 
1 .5007-01-2.2007*01 


1 

5 

If $1  i 97  I 7 3941 9! 

15  4 6.330! 

7757  7 Ml 
7737  » k#l!t 
ff  S 7 PI  J 7 lim 
?.94J»-07  7 .5?  71 -06 
1.25JI-C7  1.2577  -04 
1 .505* -07  7.9417-04 
I . 15  2* -07  1 . 7271  -04 
'•  .2007  -0  7 1.4507  -04 
6.4001-07  1.2501-04 
6.9167-0 7 7.6611-06 
5.7017-0?  5. >667  -06 
9 .690* -0  7 6.7667  -06 
1 .2621-06  l .6291  -05 
1 .6  767-06  1 . 2141-05 
1.7661-36  #.#601-06 
5 .6621-06  6. #627-06 
1.2177-25  1.5661-06 
5.0067-05  1.6961-06 


1-1. 0001*06  I I 

5 1.20)7*04  6.0001*06  1.0007*16  4.0007*0? 
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DISTRIBUTION 


t 

i 

r 


» 


i 


DEFENSE  DOCUMENTATION  CENTER 
CAMERON  STATION,  BUILDING  5 
ALEXANDRIA,  VA  22314 
ATTN  DDC-TCA  (12  COPIES) 

commander 

USA  RSCH  t STD  GP  (EUR) 

BOX  65 

FPO  NEW  YORK  09510 
ATTN  LTC  JAMES  M.  KENNEDY,  JR. 

CHIEF.  PHYSICS  t MATH  8PANCH 

COMMANDER 

US  ARMY  MATERIEL  DEVELOPMENT 
i.  READINESS  COMMAND 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 
ATTN  DRXAM-TL,  HQ  TECH  LIBRARY 

COMMANDER 

USA  ARMAMENT  COMMAND 
ROCK  ISLAND.  IL  61201 
ATTN  DRSAR-ASE,  FUZE  DIV 
ATTN  DRSAR-RDF,  SYS  DEV  DIV  - FUZES 

COMMANDER 

USA  MISSILE  6 MUNITIONS  CENTER 
< SCHOOL 
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